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ABSTRACT 

We consider the alignment torque between a spinning black hole and an accretion 
disc whose angular momenta are misaligned. This situation must hold initially in al- 
most all gas accretion events on to supermassive black holes, and may occur in binaries 
where the black hole receives a natal supernova kick. We show that the torque always 
acts to align the hole's spin with the total angular momentum without changing its 
magnitude. The torque acts dissipatively on the disc, reducing its angular momentum, 
and aligning it with the hole if and only if the angle 9 between the angular momenta 
Jd of the disc and Jh of the hole satisfies the inequality cos6> > —Jd/Uh- If this 
condition fails, which requires both 9 > tt/2 and Jd < ZJh, the disc counteraligns. 
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1 INTRODUCTION 

In a recent paper, Volonteri et al (2005; see also Madau 
2004) consider how the spins of supermassive black holes 
in galaxies change as the holes grow through both mergers 
with other holes, and gas accretion. Mergers occur at ran- 
dom angles, and when integrated over the mass distribution 
expected in hierarchical models lead neither to systematic 
spin-up nor spin-down. Gas accretion, driven for example 
by minor mergers with satellite galaxies, is likely also to oc- 
cur at random angles, and thus be initially retrograde with 
respect to the hole spin in half of all cases. However, Volon- 
teri et al. argue that gas accretion nevertheless produces sys- 
tematic spin-up, because a black hole tends to align with the 
angular momentum of an outer accretion disc on a timescale 
typically much shorter than the accretion timescale for mass 
and angular momentum (Scheuer & Feiler, 1996: hereafter 
SF96; Natarajan & Pringle, 1998). Volonteri et al. (2005) 
note that this conclusion holds only if most of the accretion 
takes place through a thin accretion disc. If instead accre- 
tion is largely via a geometrically thick disc (as happens if 
most mass accretes at super-Eddington rates) alignment oc- 
curs only on the mass accretion timescale. In this case there 
would be no net long-term spin-up, assuming successive ac- 
cretion events were randomly oriented. 

Here we address the uncertainties in our current un- 
derstanding of the evolution of warped accretion discs, and 
the resulting uncertainties in the alignment mechanisms and 
timescales for discs and black holes. 

We stress that throughout the paper we neglect the 
change of the black hole spin as it gains mass from the disc, 



i.e. we consider timescales shorter than that for increasing 
the black hole mass significantly. Thus all the torques we 
consider arise from the Lense-Thirring effect on a misaligned 
disc. (These torques are dissipative, and can cause changes 
in the local accretion rate in the disc. However these changes 
are at most by factors ~ 2.) 

We find that under some conditions counteralignment 
occurs, contrary to what is usually thought. As mentioned 
above, initial misalignment must characterise most gas ac- 
cretion events on to supermassive holes. Stellar mass black 
holes accreting from a binary companion may also be mis- 
aligned in cases where the hole received a supernova kick 
at formation. We begin by summarizing briefly what is cur- 
rently known in various cases. 



2 DYNAMICS OF ALIGNMENT 
2.1 High viscosity discs 

The dynamics of the process of alignment of black hole and 
accretion disc is not fully worked out. The best-understood 
case is when the viscosity is sufficiently high and/or the disc 
sufficiently thin that its tilt or warp diffuses through it in the 
way envisaged by Pringle (1992). This is also the most likely 
case for black-hole discs in AGN and binary systems (Wi- 
jers & Pringle, 1999; Pringle, 1999). SF96 consider the case 
where a thin, high-viscosity disc is misaligned by a small 
fixed angle 6, and linearize in 9. In this approximation the 
diffusion equation governing the time evolution of the disc 
surface density remains unchanged to first order, and the 
analysis assumes that the disc has reached a steady state. 
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In practice this requires that we consider timescales longer 
than the inflow timescale at (and within) some relevant ra- 
dius -Rwarp- SF96 (see also Rees, 1978) find that 

Rwarp ~ W p / V. (1) 

Here uj p = 2GJ/c 2 , where J = acM(GM/c 2 ) is the angular 
momentum of the hole (with — 1 < a < 1), and v is the 
kinematic viscosity in the disc. 1 

At radii R < i? W ar P the spins of the disc elements are 
aligned with the spin of the hole J^. At radii R> Rwarp, 
the spins of the disc elements make an angle 9 < 1 to Jj,. 
A spinning black hole induces Lense-Thirring precession in 
the misaligned disc elements, and the precession rate falls off 
rapidly with radius (oc R~ 3 ). Thus this induced precessional 
torque acts mainly in the region around radius i? W ar P - In 
Cartesian coordinates let define the z-axis and consider 
an elemental annulus of the disc with spin A 3d which is not 
parallel to Jj. Then each such disc annulus feels a torque 
which tries to induce precession about the z-axis. Integrat- 
ing over all these torques to get the net torque on the disc 
(and by Newton's third law the net torque on the hole) we 
conclude that the net torque can only have components in 
the [x, j/)-plane. In the analysis of SF96 they assume that 
the disc extends to infinite radius, and thus that the angu- 
lar momentum of the disc, 3d, dominates that of the hole. 
If we assume that the disc is tilted such that 3d lies in the 
(x, z)-plane, then the conclusion of SF96 is that the x and y 
components of the torque are equal in magnitude. The sign 
of the y component, which gives rise to precession, depends 
on the sign of the z component of 3d- The x component af- 
fects the degree of misalignment of the disc. It is negative 
both when the disc and hole are nearly aligned (as above) 
and also when the disc and hole are nearly counteraligned 
(i.e. when • 3d < ). This implies that the net result is 
to try to align the spins of the disc and the hole. Thus the 
prediction of the SF96 analysis appears to be that eventu- 
ally 3h and 3d should end up parallel. We see below that 
this conclusion does not hold in general. 

2.2 Low viscosity discs with small a black holes 

If the disc is thick and/or its viscosity low, so that a < H/R, 
then its tilt or warp propagates as a wave rather than diffus- 
ing (Papaloizou & Lin 1995). If such a disc is nearly aligned, 
relativistic precession effects do not align its inner regions 
with the symmetry plane of the black hole, in contrast to 
the viscous disc (a > H/R). Instead the disc tilt oscillates 
(Ivanov & Illarionov, 1997; Lubow, Ogilvie & Pringle, 2002) 
with an amplitude proportional to R 1 ^ 8 / (EH) 1 ^ 2 . In the 
inner regions of black hole accretion discs this quantity typ- 
ically increases with decreasing radius (Shakura & Sunyaev, 
1973; Collin-Souffrin & Dumont, 1990). Thus even if the 
degree of misalignment is small in the outer disc it can be 
large in the inner disc, and the angular momentum vector 
of the matter actually accreted at the horizon can make a 
large angle to the hole's spin. If the disc and hole are close 
to being counteraligned, Lubow et al. (2002) show that the 
inner disc and the hole align. 

1 SF96 assume that the effective viscosities in the disc, v\ and U2 
(Pringle 1992), are comparable. 



What interests us here are the disc torques on the hole. 
In Appendix 1 we summarise the particular analytic solution 
for the zero viscosity case presented in Lubow et al (2002) 
and give explicit expressions for the torques both when the 
disc and hole spins are nearly aligned, and when they are 
nearly counteraligned. In the absence of viscosity the torques 
result in mutual precesion of the disc and the hole. In Ap- 
pendix 2 we show how this analysis can be modified to take 
account of a small viscosity, and we investigate the modified 
torques when this has been introduced. In the nearly aligned 
case, we find that the magnitudes of the torques are sensitive 
functions of the exact disc parameters, because of the oscil- 
latory nature of the disc tilt in the inner regions. However, 
while the effect of this is to introduce great uncertainty into 
the direction of the component of the torque which gives rise 
to precession, we find that the sign of the component of the 
torque which affects the angle between the spins of the hole 
and the disc is exactly the same as in the viscous case anal- 
ysis by SF96. We note that in practice this small-viscosity 
case is most unlikely to hold for the alignment process in 
most black-hole discs in AGN or X-ray binaries (Wijers & 
Pringle, 1999; Pringle, 1999) but that warp waves may be 
important in the centres of these discs, where H/R ~ 1. 



3 GENERALIZATION 

Both sets of analysis reported above assumed that the degree 
of misalignment was small (that is, the hole and disc were 
either nearly aligned or nearly counteraligned) and that the 
disc tilt remained fixed at large radius (that is, the angular 
momentum of the disc dominates that of the hole) . We now 
consider the physics of the general case where the angle of 
misalignment is not assumed to be small, and the disc tilt 
is not assumed to be fixed. The angular momentum of the 
hole 3h is well defined. We shall denote the angular momen- 
tum of the disc as 3d, but note that this is not a well defined 
quantity. We discuss the exact meaning of Jd in this context 
below (Section 4.1). From these we construct a third vector 
representing the total angular momentum 3 t — J/i + 3d, 
which is therefore a constant vector. The torques we are 
interested in come about solely because there is a misalign- 
ment. We now define the misalignment angle 8 by 

cos6» = 3 d -3 h , (2) 

where 'hat' indicates a unit vector. We define 8 so that 
< 8 < it, with 8 = corresponding to full alignment and 
8 = 7r corresponding to full counteralignment. The degree of 
misalignment is measured by the vector A J^, and so any 
torques (which are vectors) must depend on this quantity. 
Note that this vector is zero both for 9 = and for 8 = n. 
Then in the above discussion (Section 2.1), the y axis is in 
the direction of A 3d, and the x axis is in the direction 
3 h A (3h A J<j). 

We have argued above, and indeed the analyses of both 
the high and low viscosity cases confirm, that the torque on 
the hole cannot have a component in the direction of J^. 
Thus the torque must have the form 

^ = -jr 1 [J k AJd]-^a[JfcA(JfcAJd)]. (3) 
at 

Here the first term through the quantity K\ gives the mag- 
nitude and sign of the torque which induces precession. It 
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does not lead to a change in 9. The second term describes 
the torque which changes the alignment angle 9. Both sets 
of analysis in the high and low viscosity cases show that K2 
is a positive quantity whose magnitude is dependent on the 
properties of the disc and the hole. Indeed we show below 
1100 quite generally that K2 must be positive in the pres- 
ence of dissipation. In general of course K2 is likely to be a 
function of 9 as well. 

If we take the scalar product of this equation with Jh 
we see that dj\jdt = 0, so that the magnitude of the spin 
of the hole remains constant, i.e. Jh = constant. Thus the 
tip of the Jh vector moves on a sphere. The total angular 
momentum Jt = Jh + Jd is of course a constant vector, 
representing a fixed direction in space. Using this, and the 
fact that Jh ■ dJh/dt = 0, we see that 



Ji 



Jh + Jd 



dJ h 

At 



• Jt) = J* 
Using this leads to 



AJh 
At ■ 



-(Jh-Jt)=K 2 [J 2 d jl-{J d -J h ) 
at 



! ] = A > 0. 



(4) 



(•>) 



Now since both Jh, Jt are constant, this means that 
^(cos0 h )>O (6) 

where 9h is the angle between Jh and the fixed direction 
Jt. Thus cos 9h always increases, implying that 9h always 
decreases. This means that the angular momentum vector of 
the hole always aligns with the fixed direction corresponding 
to the total angular momentum vector Jt. 

To see how Jd behaves during this process we consider 
the quantity A defined above (equation 5). We have, using 
the definition of J t and the fact that Jh = constant, that 



A = ^(J h • J t ) = ^{Jt + J h ■ J d ) = f t (J h 



Thus 

5<* 



J d )=A> 0. 



We have also that 



±J 2 
At 1 



^(J 2 h +2J h - J d + Jl) 



which implies that 



d 2 
Ai Jd 



-2J d 



dJ fe 
At 



-2A < 0. 



(7) 



(8) 



(9) 



(10) 



From this we conclude that the magnitude of the disc 
angular momentum jj decreases as Jh aligns with Jt. This 
is to be expected since, although of course the total angular 
momentum of the system (hole plus disc) is conserved, the 
alignment process requires dissipation. Since the magnitude 
of the spin of the hole remains unchanged, the dissipation 
must imply a reduction in (the magnitude of) the disc an- 
gular momentum. This justifies our statement above that 
K2 > 0, as a negative K3 would require energy fed into the 
disc rotation. Since we consider timescales short compared 
with that for accretion this is not possible. 

It is now straightforward to discover when the dissipa- 
tive torque leads to alignment or to counteralignment of Jh 
and Jd- Since the angle between Jh and Jd is 9, the cosine 
theorem gives 



cos 9 < 



(11) 

7r) occurs if and only if 
(12) 



2J/ l JdCOs(7T - 6) 

Evidently counteralignment (9 - 
Jh > Jt ■ This is equivalent to 

Jd 
2J h 

Thus counteralignment is a possible outcome and requires 
9 > tt/2, Jd < 2J h . (13) 

So why did the analysis of SF96, and also that given in 
the Appendices, imply that the disc and hole always ended 
up aligned? The reason is that both sets of calculations made 
the assumption that the outer disc was fixed, that is that 
Jd 3> Jh- In this case we see that counteralignment is for- 
bidden, and alignment must result. 

If instead Jd < 2Jh, then for Jh, Jd in random direc- 
tions, counteralignment occurs in a fraction 



Jd 
2J h 



(14) 



of cases. The disc spin is given by the relevant root of ill 1 1 . 



Jd = — Jh cos 9 + (J t 
for alignment, and 
Jd = - Jh cos9 - (J t 2 



JhSin 9) 



1/2 



Jh sin 2 1 



vl/2 



(15) 



(16) 



for counteralignment. In both cases J d decreases monotoni- 
cally in time, reaching the final values J t — Jh and Jh — Jt 
respectively. 

We can now derive the equation governing the change 
of 9. From JH) we have 



A 



Jh (Jd cos ( 



JhJd-TT (COS 0) 



A^-cos9 

Jd 



(17) 



where we have used ll 1 ti to write AJd/At — —A/J d . Collect- 
ing terms and noting that A — KiJhJd sm2 ^ we have 



4 (cos (9) 
AV ' 



KiJh sin 2 9{J d + Jh cos 9) 
and from 1151116(1 we get 

±K 2 J h sm9{Jt 



4- (cos 6») 



Jh sin 2 f 



,1/2 



(18) 



(19) 



where +/— corresponds to alignment/counteralignment re- 
spectively. Expanding these two equations about 9 = 0, n 
respectively shows that these equilibria are stable in the two 
cases. Note that there is no contradiction between the global 
alignment criterion ll I ill and the local equation l|18|l : 9 does 
not always decrease monotonically for alignment or increase 
monotonically for counteralignment. 

3.1 A geometrical picture 

Although in terms of the algebra given above the alignment 
process looks somewhat complicated, in terms of the geom- 
etry of the situation it is quite simple. 

In Figure^a) the initial vector Jd is represented by the 
line OA, and the initial vector Jh by the line OB. Then the 
total angular momentum Jt is represented by the line OC, 
where OACB forms a parallelogram. Thus throughout the 
subsequent evolution the line OC remains fixed. Since we 
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Figure 1. The evolution of hole and disc angular momenta 
Jh(OB),J £ ; (OA) under the alignment torque, viewed in the plane 
they define. This plane precesses around the fixed total angular 
momentum vector Jt (OC). (a) A case where the initial angle 6 
between Jh, 3d satisfies cos£> > — Jd/2Jh : the two angular mo- 
menta align, (b) A case where cos 9 < — Jd/2Jh'- the angular 
momenta countcralign. (c, d) Two cases where Jd 3> Jh as con- 
sidered by SF96 for which alignment always occurs. 



are just interested in the alignment process, rather than any 
precession around J t , we need only consider what happens in 
the plane defined by OACB. In this plane we have seen that 
Jh remains constant, and that the effect of the evolution is to 
align Jh with Jt. Thus, as shown in the Figure, the tip B of 
the vector OB describes the arc of a circle centred on O and 
ending up on B'. Once full alignment has occurred, the final 
vector Jh lies along OB'. Then in order that total angular 
momentum be conserved, the tip A of the vector must 
move along a corresponding arc, centred on C, and ending 
at A'. Note that as this occurs, Jd decreases monotonically. 
We see that in this case (Figure 1(a)) the final vector Jd lies 
along OA', and the disc and hole end up aligned. 

In Figure 0^b) we show exactly the same procedure, 
but with different initial values for Jd and Jh- As before 
the vector moves from the initial position OB along an 
arc centred on O to a final, aligned position, OB'. The total 
angular momentum Jt, represented by OC, remains fixed. 
The disc angular momentum vector OA moves along an arc, 
centred on C, to its final position OA'. But now, because of 



the initial values of Jd and J^, the disc and the hole end up 
counteraligned. 

In Figures^c) and^d) we show the same evolution but 
in the case Jd 2> Jh considered in the analytic calculations 
of SF96 and in the Appendix. We can see here that both disc 
and hole always end up aligned, independent of the initial 
alignment of the hole relative to the disc. 

3.2 Variation of 9 with time 

Figure |21 illustrates the evolution of misalignment angle 6 
for a particular evolution rate, namely for K2 constant in 
equation JTHJ. The general behaviour of the solutions is in- 
dependent of the detailed form of K2, provided that K2 is 
positive (i.e., dJd/dt < 0). 

The highest and lowest curves in Figure |5|show a mono- 
tonic approach to misalignment and alignment, respectively. 
The middle two curves are close to the bifurcation between 
the two end states. From equation 1181 . it follows that non- 
monotonic behaviour in time can occur for the case of a 
misaligned disc that begins an approach toward alignment, 
but in the end becomes misaligned (as seen in the sec- 
ond highest curve of Figure |2J . This situation is realised 
for an initial state having Jd < — 2 A cos f? < 2Jd- But in 
such cases, the disc never gets close to alignment; that is, 
6 never drops below n/2. A bifurcation between end states 
of alignment/misalignment occurs for cases where initially 
Jd = —2Jh cos 6. A misaligned disc that changes its sense of 
evolution (i.e., the sign of 8) from heading toward alignment 
to heading toward misalignment does so once Jd has been 
sufficiently reduced. The extent of the required reduction is 
larger for initial values of Jd/Jh closer to the bifurcation 
value — 2cos#, for which misalignment is achieved at a time 
when Jd(t) is zero. 

Figure [3] illustrates the critical case of bifurcation, ac- 
cording to the geometrical picture of the previous subsec- 
tion. Since Jh = Jt in this case, we see that as J^ aligns 
with J t , Jd approaches zero. 

Non-monotonic behaviour in the opposite sense from 
what is discussed above, i.e., starting with an evolution to- 
ward alignment and ending with misalignment, is only pos- 
sible if Jd increases in time. As described in Section 4.2, 
this can happen because of the uncertainty of what what is 
meant by Jd in a real disc. Examples of such cases are de- 
scribed in Section 4.2. These cases involve the complication 
that Jt also changes in time. 

3.3 Comparison with numerical simulations 

To date there are two numerical simulations in the litera- 
ture which have Jh > Jd/2 and thus potentially allow the 
counteralignment we predict. Nelson & Papaloizou (2000) 
perform a quasi-Newtonian 3D SPH simulation in which 
Jh >> Jd- They show that an initially prograde disc does 
align with the hole, effectively also demonstrating again 
that K-2 > in the presence of dissipation. However they 
do not consider any cases with an initially retrograde spin 
(6 > tt/2). 

Fragile & Anninos (2005) give fully relativistic 3D grid- 
based simulations, again with A >> Jd- However there 
is no explicit dissipation in their code (K2 — 0), so any 
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Figure 2. The time evolution of the disc-black hole misalignment 
angle 8 in degrees as a function of dimcnsionless time, which 
is normalized by r, the disc spin-down timescale for 8 = 90°. 
The evolution is determined by equation 1181 with an assumed 
constant value of K2 = 1/(tJ^). Initial misalignment angles 
are 9 = 135° for the upper three curves and 90° for the low- 
est curve. The initial angular momentum ratios from the high- 
est to lowest curves are Ja/Jh = 0.5,1.40,1.42 and 0.5, respec- 
tively. The curves, from the highest down, correspond to the cases 
a, b, d, and c, respectively in Fig Q For initial misalignment 
angle 135°, equation 1121 predicts that the transition between 
long-term alignment and counteralignment occurs when initially 
Jdl ' Jh = %/2 ~ 1.414, as displayed in the middle two curves that 
are on opposite sides of the transition. Notice that the second 
highest curve shows non-monotonic behaviour in time. 



A 




Figure 3. A geometrical picture of the critical case at the 
boundary between alignment and misalignment, following the no- 
tation of Figure Q In this case, we have Jt = Jh, with initially 
Jd = —Uh cos S. The end state has = Jd and Jd = 0. 

(counter)alignment can only occur on a long timescale asso- 
ciated with numerical dissipation. There is indeed a hint of 
counteralignment in the near-retrograde case they discuss. 



4 DISCUSSION 

We have considered the interaction between a misaligned 
accretion disc and a rotating black hole. We have argued 
that the torque between them must have the form JJ. The 
net result is twofold. First, there is a component (Ki) which 
causes the disc and hole to precess around the direction of 
the total angular momentum vector. The direction and rate 
of precession can depend in a complicated way on the prop- 
erties of the disc. Second, there is a torque (K2) which, since 
K2 > 0, acts to align the hole with the total angular mo- 



mentum without changing its spin rate. This torque acts 
dissipatively on the disc, and counteraligns or aligns it with 
the hole according as the conditions il l Ml hold. 

In the high viscosity case most relevant for black-hole 
discs the steady disc shape is relatively simple (SF96). In 
both the co- and counter-rotating cases the disc is flat but 
inclined to the hole at large radii, and flat but aligned with 
the hole at small radii. The change between the two (the 
warp) occurs at a radius -Rwarp given by where the rate 
at which the disc is twisted, i.e. the Lense-Thirring preces- 
sion rate ui p /R 3 , is balanced by the rate ~ v/R 2 at which 
viscous torques can propagate the twist away. Here v is the 
viscosity relevant to the process of smoothing disc warp. 
It corresponds to the viscosity V2, introduced by Pringle 
(1992), which measures the viscosity corresponding to the 
(R, z) -component of the stress tensor. 

The actual dynamics of the various alignment processes 
are likely to be complicated and need further investiga- 
tion. As we have seen in Section 3, the torque acts dissi- 
patively on the disc, reducing its angular momentum. If the 
total disc angular momentum is Jd then the disc eventually 
aligns with the hole if and only if the angle 6 between the 
angular momenta Jd of the disc and of the hole satis- 
fies cos 9 > —Jd/2Jh- If this condition fails, which requires 
6 > 7r/2, Jd < 2Jh, the disc eventually counteraligns. This 
result just follows from the physical nature of the torques, 
together with the fact that the process is dissipative, so that 
K2 > 0. Based on these simple physical ideas we were able to 
sketch the evolution of the two vectors Jh and Jd- However, 
what we are not able to do, without further consideration of 
the detailed properties of the disc in the form of the coeffi- 
cient K2, is to predict the timescale on which this happens. 



4.1 The meaning of Jd 

So far we have been deliberately vague on the precise mean- 
ing of the disc angular momentum Jd- For an accretion disc 
we may define the angular momentum vector Jd(R) of the 
material inside some radius, R. As an example we consider 
the disc model for AGN discs given by Collin-Souffrin and 
Dumont (1990). For this disc model we are interested in 
the innermost region (called Regime A, which corresponds 
to region (b) in the disc models of Shakura and Sunyaev, 
1973). In this regime if we define the radius inside which 
the angular momentum of the disc equals that of the hole 
as Rj, so that Jd(Rj) = Jh, then it is given in terms of the 
Schwarzschild radius of the hole, R„, as 

D f c \ 6/19 f T \ -6/19 

|=3.9x 10 3 ( — ) f 

R s \0.lj KO.ILeJ 

XM * 12/I9 (^k) 8/I9al ° /19 - (20) 

Here e is the efficiency of the accretion process (i.e. L = 
eMc 2 ), L is the accretion luminosity, Le is the Eddington 
limit, Mg is the mass of the black hole in units of 10 s M©, 
a is the Shakura and Sunyaev (1973) viscosity parameter, 
and a is the (dimensionless) spin of the black hole. 

The timescale on which this disc radius can communi- 
cate with the central disc regions is the viscous timescale at 
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this radius and is given by 

/ , \ 16/19 / T \ -16/19 

U( RJ ) = 1.65 x 10* (-) (_) 

xMr^r^V. (21) 

Thus on timescales longer than this we expect the ef- 
fective angular momentum of the disc to dominate that of 
the hole and therefore that on long timescales the spin of the 
hole ultimately aligns with that of the disc as in Figurcs^c) 
and (d). 

But on timescales less than this, the angular momentum 
of those parts of the disc which are able to interact with the 
hole is much less that that of the hole. On these shorter 
timescales we might expect the disc evolution to resemble 
the evolution shown in Figures ^a) and (b). Thus we have 
the apparently contradictory possible scenario in which on 
short timescales the disc tries to counteralign with the hole, 
but on long timescales t 3> t v (R.j) it ends up co-aligning 
with the hole. This means that the actual disc evolution 
depends crucially on how the warp is propagated radially 
by the disc. In other words we need to be able to predict the 
nature of the (R, z)-stress denoted by the second viscosity 

4.2 Warp propagation 

If the degree of warping is very small compared to the disc 
thickness H/R then Papaloizou & Pringle (1983) showed 
that, because of resonant effects, the warp stress is much 
larger than the usual azimuthal stress. If the warp stress 
is parametrised by 0:2 and the usual viscosity by ai, then 
they found that Q2 = l/(2ai). In this case the warp ra- 
dius can be quite small — R w /R 3 ~ 10 - 100 (Natarajan 
& Pringle, 1998). However once the warp becomes signifi- 
cant the approximations made in this analysis break down. 
One possibility then is that the resonant flows become un- 
stable (Gammie, Goodman & Ogilvie, 2000), the flow be- 
comes turbulent, and 02 is reduced significantly until per- 
haps 0:2 ~Qi. If ct\ — ct2, which is the assumption made by 
SF96, we find that 




In reality, it expected that disc viscosity is generated by 
MHD turbulence, instigated by the magneto-rotational in- 
stability (Balbus & Hawley, 1991). How this turbulence in- 
teracts with a rate of strain in the (R, ,z)-direction has yet 
to be fully determined (see, for example, Torkelsson et al., 
2000). It is evident that for finite amplitude warps and mis- 
alignments, in order to estimate the timescales and mecha- 
nisms for warp propagation it will be necessary to undertake 
numerical simulations. 

One of the first goals for such simulations will be to 
determine whether the simple picture of two viscosities is 
adequate to a first approximation (cf. Ogilvie, 2000) . Even 
in this picture it is clear that a simple a prescription is inad- 
equate. For example Larwood et al. (1996) show that when 
a disc is subject to strong forced precession (as is likely to 
occur in the inner regions of a tilted disc around a Kerr 



black hole) the disc may break in the sense that the disc tilt 
shows a sharp jump at some radius. This can only happen 
in the diffusive picture if the diffusion coefficient (1/2) is a 
function of the (gradient of the) disc tilt angle. But if this 
does happen, it enhances the possibility, discussed above, of 
the inner disc regions being able to counteralign on short 
timescales, before eventually co-aligning at later times. If 
this happened, then accretion onto the hole would act ini- 
tially (on timescales t <C t u (Rj)) to spin the hole down, 
in contradiction of the assumption made by Volonteri et al. 
(2005). 

4.3 Black holes in X-ray binaries 

Maccarone (2002) reports that in at least two of the soft 
X-ray transient (SXT) binaries (GRO J 1655-40 and SAX 
J 1819-2525) the observed relativistic jets appear not to be 
perpendicular to the orbital plane. If the jet directions are 
indicative of the direction of the spin of the hole, then the 
most likely explanation is that the misalignment occurred 
during the formation process of the black hole, and that 
subsequent evolution has not had time to bring about align- 
ment. This interpretation is interesting in that it points to 
black hole formation in a (presumably anisotropic) super- 
nova explosion. 

From 12011 we see that for a M ~ 10 M Q black hole rele- 
vant for such binary systems the radius Rj is typically much 
larger than the binary separation. Thus in these systems 
Jd *C Jh- However, the angular momentum in the binary 
orbit is much larger than that of the hole. Thus the crucial 
timescale in these systems is the timescale on which tidal ef- 
fects can transfer angular momentum from the binary orbit 
to the disc. On timescales shorter than this, the evolution 
of the disc tilt follows that shown in Figures HTa - ) and|Tfb). 
with the possibility that the disc can counteralign with the 
hole. Again, numerical simulations are required to provide 
estimates of the tidal torques for strongly misaligned discs. 

Estimates of timescales from stellar evolution theory 
can thus give lower limits to the alignment timescales in SXT 
binaries. Maccarone (2002) concludes that current theoret- 
ical estimates indicate that alignment timescales are likely 
to be at least a substantial fraction of the lifetimes of these 
systems. In any case the long quiescent intervals (10-50 yr 
or more) in SXT binaries strongly suggest that the inner 
regions of the disc are either absent or very tenuous. This 
means that virtually all of the disc mass is far outside the 
warp radius (~ a few Schwarzschild radii) and so the align- 
ment torque must be very weak. 
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Appendix 1: zero viscosity discs with small a black 
holes 

Lubow et al. (2002) consider the torque exerted be- 
tween the central black hole and a disc for the case of zero 
viscosity. They use the same linearised approximation as 
SF96, with the disc again assumed steady. At each radius 
the disc angular momentum is in the direction of the unit 
vector l(i?) ~ (l x ,l y ,V), where l x ,l y <C 1. They then use 
the complex quantity W(R,t) — l x + il y to describe the 
shape of the disc. They consider explicitly a simple exam- 
ple in which the amplitude of oscillations is independent of 
radius, and for which there is a simple analytic solution for 
a? <C r = Rj (GM/c 2 ). They take the disc thickness to vary 
as 



H h -i 
R= eT ' 

and the disc surface density to vary as 
E = E D r3-\ 

where we require h > — |. 
In this case, if we set 



24 | a | \ 1/2 r-t^ 1 / 4 ' 



(23) 



(24) 



(25) 



where a is the black hole spin parameter with positive (neg- 
ative) a corresponding to alignment (counteralignment). 
If a > 0, Lubow et al. (2002) find that 



w = w a 



cos(a;i n — x) 
COS x ln 



(26) 



where Woo gives the tilt at large radius, x- m corresponds to 
the inner boundary where the torque vanishes, i.e. dW/dr — 
0, and we need the proviso that cosxin 7^ 0. 

If a < 0, it is simple to show that the corresponding 
expression is 



w = w a 



cosh(a;i n — x) 

' COsh Xi n 



(27) 



If T x and T y are the components of the torque on the 
disc, assuming as above that the hole spin is aligned along 
the z-axis, and writing T = T x + iT y , we find that 



T = 



V6in. 1/a G 2 A/ 2 E D 
— - — a e ^ 



Wdx. 



Then in the case we are considering for a > 



T tx i I a 



,1/2 



Woo tan x in 



and for a < 



T tx 



,1/2 



Woo tannin. 



(28) 



(29) 



(30) 



In both cases the coefficient of proportionality is real and 
positive. Thus, in this case with zero disc viscosity the torque 
causes a mutual precession on the hole and the disc. But 
unlike in the high viscosity case, if a > the sign of the 
precession is determined by the details of the disc properties. 

Appendix 2: small viscosity discs with low a black 
holes 

It is not possible to introduce a small constant a vis- 
cosity as a small perturbation of the Lubow et al. analysis. 
This is because at large radii, the Lense-Thirring effect goes 
to zero, and so the small viscous perturbation does not re- 
main small. But, one can apply perturbation theory for a 
spatially varying a that takes the form 

a = a. r~ Y (31) 

for the case of small viscosity. 
If we write 

C = V(l+~a<0, 

where it is understood that we take the root with positive 
real part, then in the particular case considered above the 
solutions become for a > 



(32) 



W = Wo 



cos[((x in - x)] 



(33) 



(34) 



COs(Oin) 

and for a < 0, 

w = Wao c 0S Hax- -x)] _ 

cosh((,a;i n ) 

The warps acquire a y— component that is out of phase 
with respect to Woo as a consequence of the viscosity. This 
phase shift leads to a net alignment torque (nonzero x- 
component of torque) on the disc. In the a > case, radially 
oscillatory warped waves can occur provided that x ln > 2n. 
Such waves are possible for discs that are sufficiently thin 
{•Ja/e sufficiently large). In Figure 0] we plot W as a func- 
tion of r for two sets of disc parameters that differ in the 
value of e. Notice that oscillatory behaviour occurs for the 
thinner disc. 

The torque on the disc is given by 



(35) 



T _vW 1/2 G 2 A/ 2 S tan[C:r in ] 

This torque applies to both positive and negative values of a 
through analytic continuation. It is straightforward to show 
that the x-component (real part) of the torque is negative, 
independent of the sign of a. This result implies that the de- 
gree of misalignment decreases in time, as will be discussed 
more fully in the next subsection. 
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Figure 4. The real (top) and imaginary (bottom) parts of 
W (r) /Woo as a function of radius r (in units of GM/c 2 ) for 
two disc-black hole cases. In both cases, we adopt a = 0.3, and 
h = 1. The inner radius ri n occurs at the marginally stable orbit. 
The oscillatory (wave-like) solutions are for e = 0.01, while the 
non-oscillatory solutions (plotted with the heavier lines) occur for 
e = 0.1. The viscosity parameter a = er; n /r. 



We consider discs with inner truncation at the radii of 
the marginally stable orbits. In the prograde spin case (a > 
0), the alignment torque Re(T) can undergo large variations 
as a function of parameters a and e. But, there is a well 
defined average value. 

We determine the a-averaged value of the alignment 
torque in the case of a > with a fixed value of e. For 
a ^> e 2 (i.e., Xi n 3> 1), the torque undergoes multiple local 



peaks in value where £i n (a) 
Near such points, 



mv/2 for positive integer n. 



tannin] ~ ci 



,1/2 



717T 

— + icincto 



(36) 



where d are real constants of order unity and we have ig- 
nored variations in the inner disc radius (marginally sta- 
ble orbit) as a function of a. Since these peaks are spaced 
in a by an amount C3ne 2 , the a-averaged value of q(a) = 
Im(tan[C } Xi n (a)]) can then be expressed as 



< q(a) >= - 



dz 



+ 1 



(37) 



where z = (cia 1 / 2 /e — nir /2) / \na C2). w is the peak width 
expressed in terms of z which is inversely proportional to 
a . For small a , we can take the integral limits to infinity 
and we find that g is independent of e and a . Consequently, 
we can approximate the average alignment torque for a > 



Figure 5. The negative dimensionless alignment torque (top) 
and the dimensionless precession torque (bottom) on the disc as a 
function of a > 0. The dimensional torques are recovered by mul- 
tiplying by e Woo G 2 M 2 T, /c 2 . The disc inner radius ri n occurs at 
the marginally stable orbit for each a. Parameter h is unity. Two 
cases of disc-black hole systems are plotted, corresponding to the 
two cases in Figure 3 (but with varying a). The more strongly 
fluctuating torques are for e = 0.01, while the smoother torques 
(plotted with the heavier lines) occur for e = 0.1. The viscosity 
parameter a = eri n /r. 

by taking the average value of tan[("Xin]/C m equation 135> 
to be a constant of order unity that is independent of e and 
a a . 

For a fixed value of e such that Xi n 3> 1 (so that wave- 
like behaviour occurs), the a-averaged value of the torque 
from to a is approximately 

G 2 M 2 S 



< Re{T) >r 



,1/2 



eWo. 



(38) 



This average torque is then independent of viscosity param- 
eter a . The precessional torque is of similar order. 

In Figure we plot the dimensionless torques on the 
black hole as a function of a > 0. Notice that the disc align- 
ment torque is negative, indicating that alignment occurs. 
The precession rates undergo changes in sign as the spin rate 
a changes. 

For the retrograde spin case (a < 0) with Xi n ^> 1, the 
torque follows from equation l|35|l . 

V67T /. a a 

I i. -l 

6 



,1/2 



G 2 M 2 T, 



(39) 



3 V b / 1 1 c z 
In this case, we see that the ratio of the alignment torque to 
the precessional torque is a /6. Furthermore, the precession 
does not change direction as a function of a, as was found in 
the case for a > (see Fig^J. Consequently, for low values of 
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the turbulent viscosity parameter, the alignment time scale 
can be much longer than the precession time scale. This 
situation is unlike the case for a > 0, where both time scales 
are comparable. 

For both the prograde or retrograde cases, the torque 
on the disc is exerted where x ~ 1. The torque radius in 
units of GM/c 2 is then given by 



Consequently, for the purposes of computing torques, such 
as in equation l.'ffljl . the viscosity parameter a is related to 
q = OLojr by 




(40) 




(41) 



c 



for h 



1. 
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